Inertial mass of a superconducting vortex 
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We show that a large contribution to the inertial mass of a moving superconducting vortex comes from 
transversal displacements of the crystal lattice. The corresponding part of the mass per unit length 
of the vortex line is Mi = (m 2 c 2 /647r« 2 /iA|_) ln(Ai/^), where m e is the the bare electron mass, c is the 
speed of light, a — e 2 /hc~ 1/137 is the fine structure constant, fj, is the shear modulus of the solid, 
Xl is the London penetration length and £ is the coherence length. In conventional superconductors, 
this mass can be comparable to or even greater than the vortex core mass computed by Suhl. 
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PACS numbers: 74.20.-z, 74.25.Kc, 74.25.Qt 



In this Letter we analyze an old problem of the inertial 
mass of a moving vortex line in a type-II superconduc- 
tor. While in most applications of superconductors the 
dynamics of vortex lines is dissipative, the inertial mass 
can play an important role at slow motion of vortices and 
in problems of quantum tunneling of vortices. Besides, it 
is a fundamental problem of the theory of superconduc- 
tivity, not without a controversy, repeatedly addressed 
by a number of researchers over the last forty years. We 
will show that all previous works on this subject over- 
looked what appears to be a significant contribution to 
the vortex mass in conventional type-II superconductors: 
the torsional deformations of the crystal lattice by the 
moving vortex. The brief history of the subject is out- 
lined below. 

In 1965 Suhl [ [l| computed two contributions to the 
vortex mass: the mass due to the kinetic energy of 
the vortex core, M c = (2/7r 3 )mfci?, and the mass due 
to the electromagnetic energy of the vortex, M em = 
(ri)/^) 2 A/ c . Here m is the effective electron mass, ftp 
is the Fermi wave vector of normal electrons, rr) is the 
Debye charge-screening length, and £ is the supercon- 
ducting coherence length. The vortex core mass arises 
from the change in the effective mass of electrons by 
(A/ej?)m when one moves away from the center of the 
core, with A and cf being the superconducting gap and 
the Fermi energy respectively (see, e.g., Ref. [ 0). As to 
the electromagnetic mass, it is simply due to the static 
energy of the magnetic flux carried by the vortex (see, 
e.g., Ref. [[jj). The condition £ > r D yields M c > M em . 
Suhl's results were later reproduced by other methods 
(see, e.g., Refs. [ H, 0, 0, |l 0j). A separate question, 
not addressed in this Letter, is the dynamic mass of the 
vortex in a superclean limit which may be relevant to 
the low-temperature physics of clean high-temperature 
superconductors. This mass arises from the quantization 
of the electron states inside the vortex core, which is im- 
portant at small £. It was shown to exceed the core mass 
by a large factor (e F /A) 2 , [@,||. 



In this Letter we will demonstrate that a moving vor- 
tex produces torsional shear deformations of the crystal 
lattice that contribute to the vortex mass. Such deforma- 
tions are described by the transversal displacement field 
u(r, t) satisfying 



V • u = 



(1) 



In the long-wave limit they do not affect the density of the 
ionic lattice. Notice that the contribution to the vortex 
mass due to longitudinal elastic deformations (Vxu = 0) 
has been studied before [ [l(], [ll], [l2| . That contribution 
comes from the tiny difference between the mass den- 
sity of the normal and superconducting phase and has 
nothing to do with the effect studied by us. The latter 
effect arises from the electromagnetic coupling between 
the superconducting current and the local rotations of 
the crystal lattice. 

We first recall a well-studied problem of a uniform ro- 
tation of a superconductor at an angular velocity fl. A 
simple argument shows that the electron superfluid can- 
not remain at rest when the ion lattice is rotating [ [l3| . 
That would result in a huge current of ions relative to 
the electrons, and the corresponding huge magnetic mo- 
ment, M = — (e/4c)NR 2 n, where e < is the charge 
of electron, N is the macroscopic total number of super- 
conducting electrons, and R is the radius of the cylin- 
der. Obviously, this is not the case. As the ions begin 
to move with respect to the electron superfluid, the in- 
creasing magnetic flux due to the ionic current induces, 
through Faraday's law, the eddy electric field that drags 
superconducting electrons together with the rotating ion 
lattice. Some effect of the rotation remains, however. 
It comes from the Coriolis force on the electrons in the 
rotating frame of reference (see, e.g., Ref. [ 14 and ref- 
erences therein). Another way to look at it is that, in 
the rotating frame, the electrons feel the effective mag- 
netic field, H = (2m e c/e) ft. Consequently, due to the 
Meissner effect, the rotating cylinder develops the mag- 
netic moment 47rM = — (2m e c/e) J7, so that the total 
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field in the rotating frame, H + 47rM, is zero. This is 
a well-known effect predicted by London [ [l5|] and first 
observed by Hildebrandt [ |ll|. Our goal is to study the 
effect of local rotations, 



n(r,t) = ^Vxu , 



(2) 



caused by the transverse elastic deformations of the crys- 
tal lattice. 

A clear illustration of the problem in hand can be 
best achieved within the Lagrangian formalism at T = 0. 
Consider a simplified model via London [ ^5] , in which 
each ion supplies Z itinerant electrons, all of them su- 
perconducting at T = 0. In the long- wave limit the La- 
grangian of the electronic superfluid and the ionic lattice 
is 
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(3) 



Here n and n/Z are, correspondingly, the concentrations 
of electrons and ions; — Ze > and M, are the charge 
and the mass of the ion, respectively; Vi (r) is the veloc- 
ity of the superfluid in the laboratory coordinate frame, 
A is the vector potential, Xikim is the tensor of elastic co- 
efficients, and Uik = \(diUk + dkUi) is the strain tensor. 
The full Lagrangian is obtained by adding to Eq. the 
Lagrangian of the electromagnetic field, 



8vr 



d 3 r (E 2 - B 2 ) . 



(4) 



Here E = — A/c is the eddy electric field, satisfying 
V-E = 0, VxE = -B/c, and B = VxA is the mag- 
netic field, satisfying VxB = (47r/c)j + E/c, where j is 
the superconducting current. If j depends on time, the 
term — (Ze/c)u-A in Eq. (||) produces the force — ZeE 
in the equation of motion for the ions that results in the 
elastic strain. 

Since the electric current, 



. sc , . , 

j = — = env = en(v L - u) , 
oA 



(5) 



is determined by the motion of electrons with respect to 
the ions, it is natural to switch in Eq. (^) to the variable 
v = Vi — u . Then Eq. (H) becomes 



C = 



d 3 7 



d d r- 



- v- A + m e v-u 

c 



(Mi + Zm e )u 2 



^iklm^ik^lm 



(6) 



The important feature of this equation is the term 
J d 3 rnm e v-u describing the interaction between the 
electric current j = env and the dynamic transversal 



deformations of the ion lattice. We see that u plays the 
role similar to that of the vector potential A. Combining 
together the second and the third term in Eq. (^), one 
can introduce the effective vector potential seen by the 
electrons, 



A eff = A H u 

e 



and the corresponding effective magnetic field, 

B eff 



„ 2m.c„ 

VxA e// = B + — —ft 



(7) 



(8) 



where fi(r, t) given by Eq. (g) is the angular velocity of 
the local rotations. The second term in Eq. (||) is respon- 
sible for the London's magnetic moment of a rotating 
superconductor discussed above. 

In this Letter we want to study the effect of the 
transversal displacements produced in the ion lattice by 
the electric field of a moving fluxon. For that purpose, 
let us concentrate on the part of the total action that 
contains u, 
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dt 



d 3 r 



m e 
e 



1 



J u 
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(9) 



Here p = (n/ z)Mi+nTa e and \x are, respectively, the total 
mass density and the shear modulus of the solid which, 
for simplicity, we consider isotropic. The constant can 
be presented as p — pcf, where Ct is the speed of the 
transverse sound. Notice that the elastic energy of an 
isotropic solid also contains the term ^Ait 2 fc with A + |/i 
being the compression modulus [ ^1. In our case this 
term is zero due to the condition (|l]). According to Eq. 
(it), u satisfies 



d 2 u 

dt 2 



c 2 V 5 



m e dj 
ep dt 



(10) 



Looking at Eq. (|T^), one makes a startling observa- 
tion. Since j is proportional to the electric charge, the 
charge cancels out from Eq. (|lC)|), as well as from Eq. (jTJ), 
despite the electromagnetic nature of the forces acting 
between the electronic and ionic subsystems. Eq. (JTJj), 
therefore, must be the consequence of the momentum 
conservation rather than the result of a particular mech- 
anism of interaction between the particles of the crystal. 
This also follows from the fact that the first term in Eq. 
(0) can be written in the Coriolis form, J d 3 rh-fl, where 
L = (m e /c)rxj is the density of the angular momentum. 
The latter is true, however, only if m e and e in Eq. (||) 
are bare electron mass and charge, which was postulated 
by our simplified model. We shall now prove that in the 
low-frequency limit, Eq. (fJ^j, with bare electron mass and 
charge, is correct regardless of the model. 

We first notice that the electric current, expressed in 
terms of the phase tp of the superfluid wave function, is 



he / 2e A 

J= 2n^ n r^^ Ae// 



(11) 
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where A. e ff is given by Eq. (0). While ([□]) can be for- 
mally obtained from Eq. (|3|) in the spirit of the Ginzburg- 
Landau theory [ [l8| , it is in fact the most general form of 
the gauge invariant superconducting current, with m be- 
ing the effective electron mass and n being the concen- 
tration of superconducting electrons. A crucial observa- 
tion in explaining the value of the London's moment is 
that the electrons see the effective vector potential of Eq. 
(0), with m e and e, being bare electron mass and charge. 
The latter follows from the fact that equations (0) and 
(||) are the consequence of the Larmore theorem which 
provides the relation between $7 and B (or, equivalcntly, 
the relation between the mechanical angular momentum 
and the magnetic moment of electrons) in terms of the 
bare electron mass and charge regardless of the interac- 
tions. 

With Eq. ([ll]) in mind, let us study the Lagrangian of 
the system, C — C e i + £/, expressed in terms of coor- 
dinates of individual electrons r a , the phonon displace- 
ment field u, and the electromagnetic field A. Here £ e i 
includes all terms that depend explicitly on r Q , while C f 
describes long-wave deformations and electromagnetic 
fields. Consider now the Raus function [ pjJ] 



C = H e 



(12) 



which is the Hamiltonian for electrons and (with minus 
sign) the Lagrangian for u and A. Here p a = dC/dr a 
are canonical electron momenta. Because electrons adi- 
abatically adjust to the low-frequency long-wavelength 
displacements of the lattice and the electromagnetic field, 
one can quantize TL e i at any point in a solid as if A and 
u were constant. Since the electrons see the effective 
vector potential of Eq. (^|), the resulting free energy, F, 
must depend on u and A through the combination (Q), 
F = F[A e ff]. The gauge invariant superconducting cur- 
rent must now be obtained as [ |2(J 



J = -c 



8F[A eff ] _ e SF[A 



effl 



SA m e Su 

Consequently, the effective Lagrangian for u and A is 



(13) 
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eff 



-Reff = C f~ F l A eff] = C 
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e 



(14) 



The part of this Lagrangian that depends on u deter- 
mines the effective action for the deformations, 



S e// ful = I dt I d 3 r 



ny 
e 



ydu+-pu 2 -im\ k 



(15) 



The variation of S e f / [u] on u gives Eq. (|T(]) , which con- 
cludes the formal proof of that equation for arbitrary in- 
teractions between the particles of the solid. Eq. is 



exact as long as the lattice deformations are sufficiently 
slow to allow for local thermodynamic equilibrium in the 
electronic subsystem. 

We shall now proceed with the computation of the ef- 
fective mass of a moving vortex. According to Eq. (|Tc|), 
such a vortex produces torsional deformations of the crys- 
tal. If the magnetic flux was not quantized, the full treat- 
ment of the problem would require a self-consistent solu- 
tion of the equations of motion for the deformation field 
u and the vortex current j. The quantization of the flux, 
however, makes the vortex a robust source of deforma- 
tions in the right-hand-side of Eq. (|l0|). Let the vortex 
line be oriented along the Z-axis. Then, the current de- 
pends on coordinates and time through j(r — Vi), where 
r is the radius- vector perpendicular to Z . For the vortex 
moving at a speed V <C Ct, the terms in the action and 
in the equation of motion, that are quadratic on u, have 
a (V/ct) 2 smallness compared to other terms and can be 
neglected. In that limit, Eq. (^) and Eq. (|l5|) are equiv- 
alent. Integrating the first term in Eq. (|9|) by parts and 
omitting small terms, we get 



S= dt 



^u.(V.V)j-/i« 



(16) 



Since both V and j arc perpendicular to Z , it is clear 
that the spatial Fourier harmonics of u that couple to j 
must have k and u also perpendicular to Z. The problem 
then becomes essentially two-dimensional as, of course, 
is expected from its symmetry. In terms of the Fourier 
transforms of j(r) and u(r) one obtains 



S 



dt / dz 



d 2 k 



[-/Kfc 2 U k -u£ 



-5U(k.V)(jk< 

e 



Jk'Uk) 



(17) 



The minimization of this action on u k is equivalent to 
taking the Fourier transform of Eq. ( |l0| ) in the limit of 
small V. Substitution of the solution of that equation 
back into Eq. (17) yields 



S = dt dz 



dt 



MjV' 
2 

dz 



d 2 k (k-V) 2 



(2tt) 2 k 2 



-Jk'Jk 



(18) 



Further averaging over the angles gives for the vortex 
mass per unit length along the Z-axis: 



Mi = 



2e 2 fj, 



d\j 2 {r), 



(19) 



where j(r) is the current of a stationary vortex. 

The integral in Eq. @ can be expressed in terms of 
the unit-length energy of the vortex line, E v , [ ^8). In 
the limit of k = Xl/S, ^> 1, 



2tt\ 2 l 



d 2 rf(r) 



WW 



In k 



(20) 
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where $o = ch/2e is the flux quantum and \l is the 
London penetration length. This gives for the vortex 
mass per unit length 



Mi = 



n {fin) 2 

In k = 

4 /i 



— I In k , 
m 



(21) 



where we introduced the fine structure constant, a = 
e 2 /Tic~ 1/137, and used the relation \\ = mc 2 /47rne 2 to 
obtain the second of Eq. . 

Three observations are in order. Firstly, the non-zero 
value of Mi is due to the finite rigidity of the crystal with 
respect to the shear stress. In an absolutely rigid crystal, 
|U = pc 2 = oo and Mi would be zero. Secondly, M; scales 
as the square of the superfluid density, which is a rather 
unique feature provided by the specific mass-generating 
mechanism studied in this Letter. Consequently, Mi 
should go to zero as (T c — T) 2 when temperature ap- 
proaches the critical temperature T c . Thirdly, Mi does 
not fall into the category of vortex masses [ ^ that satisfy 
Mj ~ Ei/c 2 , where Ei is a contribution of some mecha- 
nism to the energy of a static vortex and Cj is the velocity 
of propagation of the distortion in question. This is be- 
cause no static energy is associated with the mechanism 
that generates Mj. 

Eq. (|2l]) shows that the vortex mass studied in this 
Letter is important in metals with high concentration of 
superconducting electrons. In copper oxides, because of 
low n, and also due to strong dynamic effects associated 
with the quantization of the electron levels in the vortex 
core [ fjj ||, g], this mass should not play any significant 
role. In good metals, however, Mi can be the main contri- 
bution to the inertial mass of the vortex line. Taking for 
estimates, n~ 10 23 cm~ 3 , p, ~ 10 11 g/cm 3 , and m~m e , 
one gets from the second of Eq. (Ell), Mi ~ 10~ 19 g/cm. 
Thus, Mi can easily exceed the inertial mass of the vortex 
core computed by Suhl, M c = (2/n 3 )mk F [ § § §. It is 
interesting to note, that, within the model of an isotropic 
good metal, these two masses do not have any parameter 
smallness with respect to each other. Indeed, taking at 



T = 0, 



k 3 F /3n 2 



Mi 



pc 2 t , 



Min/Z, we get 



(-i 



M~ = T UJ W 1 ln "' 



(22) 



where we have used the plasma approximation for the 
speed of the longitudinal sound, q = (Zva/'Z>Mi) 1 / 2 vp > 
•\/2ct [ |l7j] . However, when using the parameters of real 
materials, Mi seems always to exceed M c by one- two or- 
ders of magnitude. 

Finally, we would like to comment on the contribution 
of the interaction term in equations (]|),(fL5|) to the viscos- 
ity for the vortex motion. Conservation of linear momen- 
tum prohibits the radiation of transversal phonons by a 
vortex moving at V < Ct- At V > c t , however, and 
especially close to c t , the viscosity for the vortex motion 
may, in principle, be dominated by the radiation of the 



transversal sound rather than by the conventional mech- 
anism due to the finite normal-state resistivity [ |2l[] . The 
full solution of this problem will be presented elsewhere. 

In summary, we have computed the contribution to the 
inertial mass of a moving superconducting vortex, which 
is coming from the torsional deformations of the crystal 
lattice. Rigorous solution of this problem for an isotropic 
solid has been obtained. This contribution to the vortex 
mass can dominate over all other contributions in metals 
with high concentration of superconducting electrons. 
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